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We perform detailed magnetotransport studies on two-dimensional electron gases (2DEGs) formed 
in undoped Si/SiGe heterostructures in order to identify the electron mobility limiting mechanisms 
in this increasingly important materials system. By analyzing data from 26 wafers with different 
heterostructure growth profiles we observe a strong correlation between the background oxygen 
concentration in the Si quantum well and the maximum mobility. The highest quality wafer supports 
a 2DEG with mobility fi = 160,000 cm^/Vs at a density n = 2.17 xlO^^/cm^ and exhibits a metal- 
to-insulator transition at a critical density ric = 0.46 x 10^^/cm^. We extract a valley splitting Av 
~ 150 /xeV at a magnetic field B = 1.8 T. These results provide evidence that undoped Si/SiGe 
heterostructures are suitable for the fabrication of few-electron quantum dots. 

PACS numbers: 73.21.Fg, 73.21.La, 85.30.De 


INTRODUCTION 

The development of silicon quantum devices has gained 
considerable momentum due to reports of quantum co¬ 
herence times (T 2 ) as long as 39 minutes [T] . Its naturally 
abundant isotope, ^^Si, carries zero nuclear spin, reduc¬ 
ing hyperfine-induced dephasing due to fluctuations of 
the nuclear spin bath US]. Its small spin-orbit cou¬ 
pling is also beneficial for spin qubits Ei- Follow¬ 
ing work in GaAs quantum dots, early experimental ef¬ 
forts were made towards fabricating Si quantum dots in 
modulation-doped Si/SiGe heterostructures, where the 
n-type dopant layer is separated from the Si quantum 
well (QW) by a setback distance ranging from 5 to 20 
nm CHini- Doped devices encountered challenges when 
operating in the few-electron regime, displayed hysteresis 
in gate voltage sweeps [TO] , and sometimes suffered from 
leakage between the 2DEG and depletion gates 13 Ej. 

It is now widely accepted that the elimination of the n- 
type dopant layer decreases the Coulomb disorder in the 
QW, and reduces hysteresis and gate leakage [TTl IT2] . 
Recent experiments focusing on quantum dots made in 
undoped Si/SiGe QWs pTHT4] have consistently reached 
the single-electron regime and demonstrated inhomoge¬ 
neous spin dephasing times T 2 = 360 ns in naturally 
abundant Si, a substantial increase compared to GaAs 
spin qubits EEj. Further improvement of the Si/SiGe 
QW system may be feasible if the remaining mobility 
limiting mechanisms are clearly identified [IMT]. 

The dominant scattering sources can be identified from 
measurements of the carrier mobility /i as a function 
of 2DEG charge density n, as well as measurements of 
the quantum lifetime Tq Haim. For example, scatter¬ 
ing from remote impurities [15] is predicted to result in 
a power-law dependence yu oc Such experiments 

have been extensively performed for GaAs/AlGaAs 
heterostructures GaN/AlGaN heterostructures 

[251 [26], Si MOSFETs [27] and doped Si/SiGe het¬ 
erostructures [28H3Q] . The conclusions reached in doped 


Si/SiGe heterostructures are not directly applicable to 
undoped structures, and similar measurements on un¬ 
doped Si/SiGe heterostructures are scarce. One recent 
work reports a record-high mobility of 2 x 10^ cm^/Vs 
at a temperature T = 0.3 K and n = 1.4 x 10^^/cm^ in 
an undoped Si/SiGe QW and identifies background im¬ 
purity charges and interface roughness at the QW as the 
main mobility-limiting mechanisms [31]. However, the 
526 nm thick spacer layer used in this work is much too 
large to allow the tight electrostatic confinement that is 
needed for few-electron quantum dot devices [2] . Li al. 
investigated a heterostructure with a spacer layer thick¬ 
ness of 60 nm and deduced that remote impurity charges 
at the Al 203 /Si interface limit the mobility [30] . 

To thoroughly investigate the mobility limiting mecha¬ 
nisms in undoped Si/SiGe QWs, we report a series of sys¬ 
tematic magnetotransport measurements in the density 
range 0.5 - 8.0 x 10^^/cm^ at temperatures from 0.35 - 
4.2 K. By examining 26 different heterostructure growth 
profiles, we identify a strong correlation between back¬ 
ground oxygen concentration in the QW and maximum 
mobility. These results indicate that significant enhance¬ 
ments in Si/SiGe mobility might be obtained through 
more careful control of background contamination dur¬ 
ing heterostructure growth. 


SILICON GERMANIUM HETEROSTRUCTURES 

The samples examined here were grown at Lawrence 
Semiconductor Research Laboratory using chemical va¬ 
por deposition. 26 Si/SiGe heterostructures are inves¬ 
tigated in order to distinguish mobility limiting mech¬ 
anisms that are related to the growth profile [Fig. [^a)] 
from those that are related to background impurities. Re¬ 
laxed buffers of Sii_a,Gea, are first grown on Si substrates, 
varying x from 0 to 0.3 over a thickness of 3 /im. A 1 /im 
thick layer of SIo.tGco.s is grown on the virtual substrate 
before it is polished. The wafers are completed by grow- 
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FIG. 1: (a) Heterostmcture growth profile, (b) Top view of 
the Hall bar device. The gold region marks the area covered 
by the top gate. The other six square pads are implanted 
regions which form Ohmic contacts to the electron gas. The 
Hall bar dimensions are W = 170 jum and L = 375 jum. (c) A 
typical “turn-on” curve of the device at T = 0.35 K showing 
the source-drain current Isd as a function of Vg • The thresh¬ 
old voltage for non-zero /sd is Vr = 0.45 V. Inset shows n for 
Vg > Vt, before saturation. Dimensions and data shown in 
(a) and (c) are for Wafer No. 16. 


ing a 225 nm thick Sio.rGeo.s layer, followed by a Si QW 
which is strained through the Si/Sio.yGeo.s lattice mis¬ 
match, a Sio.yGeo.s spacer layer and a protective Si cap. 
We investigate heterostructures with Si cap thicknesses 
of 2 nm and 4 nm, Sio.yGeo.s spacer layer thicknesses h = 
20 nm, 30 nm, 40 nm and 50 nm, and Si QW thicknesses 
of 5 nm, 8 nm and 11 nm. 

Hall bars are fabricated on each of the 26 wafers, with 
the geometry shown in Fig. Bb). We first use atomic 
layer deposition to grow an AI2O3 gate dielectric on top 
of the Si cap. We then evaporate Cr/Au on top of the 
AI 2 O 3 to form a top gate. A positive dc bias is applied to 
the top gate to accumulate electrons in the QW and a 0.1 
mV, 17 Hz ac voltage excitation is applied between the 
ohmic contacts marked with S and D in Fig. [^b). The 
longitudinal voltage, I4x, and the Hall voltage, Vcy, are 
simultaneously measured as a function of magnetic field 
B using standard ac lock-in techniques. The source-drain 
current, /sd, is also measured using a current-to-voltage 
pre-amplifier and a lock-in amplifier. Employing a fixed 
voltage source in place of a fixed current source prevents a 
destructively large voltage from being applied across the 
sample at low electron densities where the longitudinal 
resistance is large. The 2D longitudinal resistivity, pxx = 
{V^^/Isd){W/L), and Hall resistivity, p^y = (V^y/IsD) 
are calculated from the measured voltages and currents. 
We have confirmed that the voltage bias is not heating 
the sample by comparing measurements of pxx and pxy at 


different excitation voltages. Density, n, and mobility, /i, 
of carrier electrons are calculated according to the Hall 
formulas n = B/{ep^y) and p = (l/5)(pxy/Pxx)- 

Figure [^c) displays a typical “turn-on” curve of the 
Hall bar devices. Zero current flow is observed below a 
threshold top gate voltage Vt = 0.45 V. For Vg > Vr, 
current starts to flow and we observe a linear increase 
in n with a slope of dn/dVc = 3.96 x 10^^/cm^/V. At 
Vg > 2.5 V, the electron density is fixed at a constant 
value of 8.0 x 10^^/cm^. Such behavior can be under¬ 
stood as follows: We define AEc = ^(Vg — Vt) to be the 
difference between the chemical potential of the 2DEG 
electrons and the conduction band minimum of 2DEG 
electrons, where e = 1.60 x 10“^^ G is the magnitude 
of the electron charge. Eor Vg < br, AEc < 0 and the 
chemical potential of the 2DEG electrons lies below the 
conduction band minimum. As a result, the QW is de¬ 
pleted of carriers. When Vg > Er, AEc > 0 and the con¬ 
duction band becomes populated. The Hall bar behaves 
like a parallel plate capacitor at these voltages. Using 
relative permittivities of = 9 for AI2O3 and = 13.1 
for Sio.7Geo.3, we calculate dn/dVc = 4.00x 10^^/cm^/V, 
which is within 1 % of the experimental value. At even 
higher values of Vg (data not shown), electrons start to 
accumulate at the Al203/Si cap interface, screening the 
QW from any further increase in Vg- This causes a sat¬ 
uration of the electron density at a constant value of 8.0 
X 10^^/cm^ for Vg > 2.5 V. Lu et al. have investigated 
such saturation behavior in detail [32] . 

CHARACTERIZATION AT T = 4.2 K 

Hall bars are first measured at T = 4.2 K and B = 
0.1 T, below the onset of Shubnikov-de Haas (SdH) os¬ 
cillations. Eigure 2(a) shows the spacer layer thickness 
h for each of the 26 wafers, along with the Si cap thick¬ 
ness and QW width. Recent studies of undoped Si/SiGe 
structures have shown that remote impurity scattering 
typically dominates in the low electron density regime, 
whereas both remote impurities and interface roughness 
dominate at higher electron densities [30l [31] . It is there¬ 
fore helpful to examine electron mobilities at both den¬ 
sity regimes. Corresponding electron mobilities are plot¬ 
ted in Eig. 2(b) for a high electron density nn = 7.0 x 
10^^/cm^ and in Eig. 2(c) for a low electron density np 
= 2.1 X 10^^/cm^. Surprisingly, both the low and high 
density mobilities show a nearly monotonic increase with 
wafer number, despite the large variation in heterostruc¬ 
ture parameters throughout this series of wafers. On top 
of this trend, abrupt dips in the mobility are observed at 
Wafer No. 20 and 26. 

Secondary ion mass spectrometry (SIMS) analysis was 
performed on each wafer to better understand the in¬ 
crease in mobility as a function of wafer number. These 
data sets are included in the supplemental material [33] . 
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FIG. 2: (a) SiGe spacer layer thickness, h, for wafers 1-26. 
Wafers with solid (hollow) symbols have a 2 nm (4 nm) thick 
Si cap. QW thicknesses are represented by symbol shapes. 
Gircles: 5 nm, rectangles: 8 nm, triangles: 11 nm. (b) /xh is 
the T = 4.2 K mobility at nn = 7 x 10^^/cm^. (c) /xl is the 
T = 4.2 K mobility at txl = 2.1 x 10^^/cm^. (d) No is the 
concentration of oxygen atoms inside the QW, obtained from 
SIMS data. SIMS is limited to measuring oxygen concentra¬ 
tions above 1 x 10^^/cm^. 


The background oxygen concentration is peaked at the 
surface of the wafer due to surface contamination and 
post-growth formation of native oxides. In addition, 
there is an oxygen peak near the Si QW, presumably due 
to the switching of the mass flow controllers in the CVD 
reactor. In Fig. [^d) we plot the concentration of oxygen 
atoms at the QW, for each wafer. For wafers 1 to 11, 
No decreases from 8.7 x 10^^/cm^ to the SIMS detec¬ 
tion threshold of 1 x 10^^/cm^. The decrease in oxygen 
concentration is correlated with the increase in mobility 
observed in Fig. 2(b-c). Wafer No. 26, which marks the 
beginning of a second cassette of wafers, shows an abrupt 
increase in Nq^ which is also correlated with a drop in the 
mobility. The combination of mobility and SIMS data 
suggest that oxygen contamination is a mobility limit¬ 
ing factor in these undoped Si/SiGe heterostructures. A 
similar correlation has been observed in undoped Si/SiGe 
heterostructures grown by molecular beam epitaxy [34]. 

In addition to the correlation between Nq and /x, the 
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FIG. 3: Wafer No. 16. pxx (red) and pxy (black) as a function 
of B, with Vg = 0.75 V and T = 0.35 K. n and fi are extracted 
from the resistivity values at B = 0.1 T. Shubnikov-de Haas 
(SdH) oscillations are visible at intermediate fields 0.3 T < 
H < 2 T. We observe clear quantum Hall plateaus in pxy at 
integer filling factors zx for B > 1.5 T. 

data show that the heterostructure growth profile im¬ 
pacts the mobility of samples later in the growth series. 
As h is increased from 40 to 50 nm for Wafers No. 15 
and 16, we observe a corresponding increase in /xh and 
/xp. For wafers 16-19, h is constant and both /xh and /xp 
show very little variation. At Wafer No. 20, h undergoes 
a large decrease from 50 to 20 nm, which is correlated 
with a large drop in mobility. For wafers 20-25, h in¬ 
creases from 20 to 50 nm and we see that the mobilities 
also recover to the values obtained from Wafer No. 19. 
It is also notable that the correlation between h and /xp 
is stronger than that between h and /xh, based on the 
relative sizes of the dips at Wafer No. 20. In contrast, 
for lower wafer numbers, the correlations between h and 
mobility are weaker, suggesting that No is the dominant 
mobility-limiting mechanism in these wafers. To obtain 
a better quantitative understanding of these correlations, 
we perform detailed measurements on Wafer No. 5 and 
16 at T = 0.35 K in order to contrast the properties of a 
low and high mobility wafer. 

HIGH MOBILITY SAMPLE 

Based on its high 4.2 K mobility, a Hall bar from Wafer 
No. 16 was cooled down in a ^He cryostat for further 
study. The oxygen content at the QW is No = 0.5 x 
10^^/cm^ [33]. Figure]^ shows characteristic plots of pxx 
and pxy as functions of 5 up to 8 T, with n = 2.17 x 
10^^/cm^. From the low field magnetotransport data we 
extract /x = 1.62 x 10^ cm^/Vs. We observe quantum Hall 
plateaus in pxy at consecutive integer filling factors v for 
B > 1.5 T, which indicates that both spin and valley 
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FIG. 4: Wafer No. 16. {i as a function of n for five different 
temperatures. The dashed lines show the predicted slopes 
for remote ionized impurity scattering, which obeys a power 
law dependence /i oc n^'^, and scattering due to impurities 
in the QW, which scales as // oc on average [15]. The 
charged impurity densities used to produce the two dashed 
lines are 3.7 xlO^^/cm^ at a distance of 50 nm from the QW 
center (where the A^Os/Si interface is) for remote impurity 
scattering, and 3.4 xlO^/cm^ in the QW center. 

degeneracies are lifted. In addition, pxx displays clear 
zeros, ruling out parallel conduction paths. For z/ > 6, 
plateaus in pxy are no longer visible, although oscillations 
in pxx are visible up to = 24. 

For a single dominant scattering mechanism, the elec¬ 
tron mobility is expected to scale as a power law of 
density m, /i oc with a scattering-mechanism- 

dependent exponent a. Figure]^ shows /i as a function of 
n at five temperatures ranging from 0.35 to 4 K. At T = 
0.35 K, /i(n) is not well described by a single exponent, 
an often observed feature in 2DEG systems [20]. Our 
data differ from previous work m, where an exponent of 
a = 1.7 is observed in the density range ofn = 0.6xl0^^ 
/cm^ to n = 4.5 x 10^^ /cm^. For n < 1 x 10^^/cm^, 
the data roughly follow the /i oc scaling, which is 
consistent with scattering due to remote charged impuri¬ 
ties US]. At higher n, /i increases at a much slower rate 
and displays signs of saturation when n > 5 x 10^^/cm^. 
The high density saturation likely arises from impurity 
charges located very near or inside the QW, which lead to 
values of /i that are only weakly dependent on n US]- It is 
notable that the mobility curves are temperature depen¬ 
dent at low densities, but all saturate to nearly the same 
high density value of 250,000 cm^/Vs. Another feature of 
the higher temperature data is that the density depen¬ 
dence of /i(n) becomes stronger, though the curvature 
persists up to 4 K. At T = 4 K, the data approximately 
follow a /i oc trend for n < 3 xlO^^/cm^. 

To further probe the scattering mechanisms that limit 
the mobility of Wafer No. 16, we measure low-field SdH 
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FIG. 5: Wafer No. 16. (a) Apxx = Pxx — pb as a function of 
1/B for three electron densities at T = 0.4 K, showing clear 
SdH oscillations. Traces have been offset by 150 Q for clarity, 
(b) Dingle plots for the data in (a). Values of Apxx used 
in this plot are the average of the maximum and minimum 
of each period of the SdH oscillations shown in (a). The 
SdH oscillations decay more slowly at higher densities, (c) 
Quantum lifetime Tq and transport lifetime rt as functions of 
n at T = 0.4 K. (d) Dingle ratio rt/xq as a function of n at 
T = 0.4 K. The dashed line is a guide to the eye showing a 
linear trend in the data. 


oscillations in the longitudinal resistivity, pxx- To facili¬ 
tate the extraction of quantum lifetimes, we subtract the 
slowly varying background from pxx as outlined by Co¬ 
leridge et al. [22] yielding Apxx = Pxx — Ph- Here pb is 
a slowly varying, polynomial background that has no de¬ 
tectable oscillatory component in the field range studied, 
indicating the lack of any appreciable parallel conduction 
path or inter-subband scattering in our sample [22] [23] . 
Apxx is plotted against 1/B for three densities in Fig. 
I^a). Clear periodic oscillations are observed, with a pe¬ 
riodicity of 4 in z^. This is consistent with the 2-fold spin 
degeneracy and 2-fold valley degeneracy at low fields. At 
higher fields B > 0.7 T, splitting of the peak in each 
period of the SdH oscillation becomes visible, which is 
a consequence of the increased Zeeman splitting which 
breaks the spin degeneracy of each Landau level. This 
splitting is examined in detail in Section VHL We ex¬ 
tract the amplitude of the oscillations in Apxx at each 
period in 1/H using the method of linear interpolation 
outlined by Padmanabhan et al. [35]. The decay in the 
resulting Apxx amplitudes is fit according to [23] : 

Apxx = 4poX(T) exp(-7r/u;crq) 


( 1 ) 
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where Tq is the quantum lifetime of the elec¬ 
trons, po is the zero-field resistivity, X{T) = 
(27r^/cB^/^<^c)/sinh(27r^/cBT/fico’c) is the temperature¬ 
damping factor, cJc = eBIvri" is the cyclotron frequency, 
and /cb is the Boltzmann’s constant. We use a constant 
effective mass m* = 0.2me, where rUg is the free elec¬ 
tron mass, for all fits [36]. Figure [^b) shows the results 
of such fits, known as Dingle plots. The slopes of the 
Dingle plots [23j are inversely proportional to the quan¬ 
tum lifetime Tq. We observe an increasing slope at lower 
electron densities, which suggests shorter quantum life¬ 
times at lower densities. We note that the Dingle plots 
show linear trends, suggesting the absence of any ap¬ 
preciable density inhomogeneities in our samples which 
would otherwise introduce uncertainties into estimates of 



In Fig. I^c), we compare the transport lifetime, rt, 
and the quantum lifetime, Tq, across the electron density 
range n = 1.8-6.8 x 10^^/cm^. Values of rt are obtained 
m from the mobility data in Fig. via rt = [inrC je and 
values of Tq are obtained from analysis of the low field 
SdH oscillations. Both lifetimes show similar dependen¬ 
cies on n and Dingle ratios, defined as range from 

4 to 7 as shown in Fig. [^d). In comparison with pre¬ 
vious work on GaAs/AlGaAs [23], GaN/AlGaN [25] and 
modulation-doped SiGe [28] where the Dingle ratios typ¬ 
ically range from 10 to 20, the Dingle ratios measured for 
this undoped sample are sizably smaller. The relatively 
small Dingle ratio indicates that large angle scattering 
plays a more dominant role in this sample than in these 
traditional systems, a situation which would arise when 
the distribution of impurities is more concentrated to¬ 
wards the location of the 2DEG m- Such a distribution 
is contrary to what is expected in an undoped system 
where charged impurities are thought to reside mostly 
in the Al 203 /Si interface [30], ^50 nm away from the 
2DEG in this sample. Our interpretation of the possible 
cause for such distribution is the peak in oxygen impu¬ 
rities at the 2DEG location [33]. Ionization of a small 
fraction of these oxygen atoms would lead to a sizable 
amount of impurity charges inside the QW, which con¬ 
tribute to large angle scattering with a Dingle ratio near 
unity. The decreasing trend of the Dingle ratio at higher 
densities also differs from theoretical calculations based 
on a single dopant sheet m and previous work on a Al- 
GaN/GaN system [25]. This deviation is interpreted to 
be due to the mitigated contribution of remote impurity 
scattering to the overall momentum scattering rate at 
higher densities, since the scattering rate oc 
for remote impurities but (x for impurities in¬ 

side the QW. At higher densities, scattering from impuri¬ 
ties inside the QW becomes more dominant than remote 
impurities, which reduces the overall Dingle ratio. 

The quantum lifetime Tq for this high mobility sam¬ 
ple is also measured as a function of T at n = 6.7 
xlO^^/cm^. Eigure|^a-b) displays low field SdH oscil¬ 



T(K) T(K) 


FIG. 6: Wafer No. 16. (a) Apxx as a function of 1/B for n 
= 6.7 xlO^^/cm^ at four different temperatures, (b) Dingle 
plots for the data shown in (a). The SdH oscillations decay 
more rapidly at higher temperatures, (c) Tq and n as a func¬ 
tion of T at n = 6.7 x 10^^/cm^. (d) Dingle ratio as a function 
of T. The dashed line is a guide to the eye. 

lations and the associated Dingle plots. The extracted 
values of Tq are plotted alongside rt, obtained by measur¬ 
ing the temperature-dependent mobility at this density 
[Fig-i c)]. The transport lifetime is relatively insensi¬ 
tive to temperature, while the quantum lifetime varies 
by nearly a factor of 4 from T = 0.4Ktol.5K. The re¬ 
sulting Dingle ratio, plotted in Fig.j^d), increases almost 
linearly with temperature from 4 to 11. These data are 
in contrast with a single-particle description of electron 
scattering in 2DEGs [HIIHIISH], where the temperature 
dependence of both lifetimes is expected to be weak in 
the T ^ Tp regime, where Tp is the Fermi tempera¬ 
ture (approximately 47 K at this density). Arapov et al. 
[39] have recently measured an InGaA/GaAs double QW 
structure and report a similar, strong temperature depen¬ 
dence of Tq at T ^ Tp, which the authors attribute to 
electron-electron interactions. A similar mechanism may 
explain the trends observed in this work. 


LOW MOBILITY SAMPLE 

We next examine data from Wafer No. 5, which has a 
much lower maximum mobility of /i = 7.5 x 10^ cm^/Vs 
at 4.2 K. Wafer No. 5 has a 2 nm thick Si cap, a h = 40 nm 
thick SiGe spacer layer, and a 5 nm wide Si QW. SIMS 
analysis shows a similar distribution of oxygen inside the 
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is likely dominated by impurity charges situated inside 
the QW, which is consistent with the higher oxygen con¬ 
tent observed in the SIMS data [33] . 

Figure [^b) shows Tq and Tt for the low mobility sam¬ 
ple at five different densities and with T = 0.4 K. Both 
lifetimes are shorter at lower electron densities, similar 
to the high mobility sample. Interestingly, despite the 
factor of ^3 difference in rt between the two samples, 
the values of Tq are very similar. This observation agrees 
well with recent theoretical results by Das Sarma et al 
[40] . who considered a two-impurity model and showed 
that increasing (decreasing) rt does not necessarily lead 
to increasing (decreasing) Tq when there is more than 
one scattering mechanism. We also plot the density- 
dependent Dingle ratio rt/rq of this sample in Fig. [^c). 
Tt/rq ranges from 1.3 at high density to 2.3 at low den¬ 
sity, significantly smaller than the high mobility sample. 
We were only able to measure Tq down to a density of 
2.85 xlO^^/cm^ in this sample due to its relatively lower 
quality which makes clear Dingle plots difficult to obtain 
at lower densities. The smaller Dingle ratio observed in 
this sample implies that the underlying scattering events 
are even larger in angle compared to the high mobility 
sample, consistent with scattering from QW impurities. 


ESTIMATE OF DEFECT DENSITIES 


FIG. 7: Wafer No. 5 (low mobility sample), (a) /i(n) at 
five different temperatures. The dashed line shows scatter¬ 
ing due to impurities in the QW (/x oc on average). The 
charged impurity density used to produce the dashed line is 
1.3 xl0^°/cm^ at the QW center, (b) Tq and rt plotted as 
a function of n, extracted from low field SdH oscillations. 
T = 0.4 K in this plot, (c) Dingle ratios rt/xq obtained from 
the data in (b). 


SiGe spacer as the high mobility sample [33]. However, 
the oxygen content in the QW is peaked at = 2.5 
X 10^^/cm^, which is five times higher than the high 
mobility sample. 

Figure shows the results of magnetotransport mea¬ 
surements on this low mobility sample. Reliable mea¬ 
surements of /i could only be performed on this sample 
for n > 1.7 xlO^^/cm^. Below this density the sam¬ 
ple is in an apparent insulating state. /i(n) is plot¬ 
ted in Fig. I^a) and increases with n, although with a 
weaker dependence than the high mobility sample. /i(n) 
is also temperature dependent, and more strongly scales 
with n at higher temperatures. At T = 0.35 K, /i(n) 
is nearly density-independent. As T increases, /i(n) be¬ 
comes more density-dependent and eventually reaches an 
approximate scaling of /i oc at T = 4 K. The smaller 
power-law exponent for this sample suggests that remote 
impurity scattering plays a less significant role compared 
to the high mobility sample. Instead, electron scattering 


Monroe et al have carefully analyzed seven scattering 
mechanisms that are potentially relevant to the Si/SiGe 
materials system [15]. Among these mechanisms, alloy 
scattering, scattering due to strain modulation, scatter¬ 
ing due to vicinal surfaces and scattering from threading 
dislocations are estimated to limit mobilities to above 
10^ cm^/Vs, two orders of magnitude higher than the 
mobilities measured in our samples. Interface roughness 
has been reported to be an important factor in a ^ 500 
nm deep Si/SiGe QW structure [31], but is expected to 
lead to a mobility that decreases with increasing density 
(a trend that is not observed in our data). We therefore 
limit our analysis to the two remaining scattering mech¬ 
anisms: remote impurity scattering and scattering from 
background charges. 

In this section we first review the theory of Monroe 
et al na. We then compare the measured mobility and 
Dingle ratio with predictions from this theory, allowing us 
to estimate charged defect densities in the high and low 
mobility samples. We limit our analysis to the T = 0.35 
K data, where T ^ Tp is satisfied throughout the density 
range studied in Fig.j^and Fig.j^a), and thermal effects 
are negligible. 

Remote impurity scattering is often identified as the 
dominant mobility-limiting mechanism in doped Si/SiGe 
heterostructures, where the dopant atoms contribute to 
the formation of a sheet of disordered charges located 
a setback distance from the 2DEG [SsHSO]- These 
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disordered charges result in potential fluctuations in the 
QW, contributing to electronic scattering and reduced 
mobility. The functional form that relates mobility to 
the electron density and 2D density of disordered charges 
Til was derived by Monroe et al.\ 

II « ( 2 ) 

where = 2 {qs = 2) accounts for the valley (spin) de¬ 
generacy. Remote impurity scattering results in a rela¬ 
tively strong density dependence fi{n) oc which is at¬ 
tributed to the fact that at higher n, the Fermi wavevec- 
tor k-p increases and electrons are scattered through 
smaller angles by the potential fluctuations of the same 
remote impurities. Monroe et al. also make a prediction 
for the density dependence of the Dingle ratio for remote 
impurity scattering: 

n/Tq « {lQTr/g^gs)zQn. (3) 

Thus, theory predicts n/rq oc n and n/rq oc Zq respec¬ 
tively. Monroe et al obtained the closed-form expres¬ 
sions of Eqns. (2) and (3) based on the assumption that 
^ 1- Since ~ 5.6 x 10^ /cm at n = 1 x 10^^ /cm^, 
this assumption translates to a setback distance zq > 18 


For the case of scattering due to more arbitrarily dis¬ 
tributed background charges where k-pZo > 1 does not 
necessarily hold, Monroe et al. derived a general expres¬ 
sion for Tt (related to g via rt = 


^ ^ S{q) 

kp Jo ^ ^1- q^/Akl ‘i.kl 


(4) 


where kp = A.'kti jg^g^. The integration variable q phys- 

ically represents the magnitude of the change in wavevec- 

tor for a given scattering event. The expression for Tq is 

2 

similar but without the angle-weighing factor 
-1 1 , S{q) 

S{q) is the power spectral density of impurity charges: 

where Qs = g2BI is the Thomas-Fermi screen¬ 
ing wavevector, = 12 is the relative permittivity of Si 
and ^ 2 D = (^v^s/27r)(m*//i^) is the 2D density of states. 
A^i(^) is the 3D density of impurity charges, where 2 : 
is measured relative to the center of the 2DEG wave- 
function ( 2 ; = 0). Based on this expression, the special 
case of a uniform background charge Ni{z) = A^b gives 
g = {gy^‘^gs^^/4:7r^^‘^){en^^‘^/hN]j). The Dingle ratio for 
a uniform background charge density was not explicitly 
derived by Monroe et al, but is expected to be large. 


Considering the SIMS analysis, which shows the pres¬ 
ence of oxygen in the QW, we also analyze scattering 
from a 2D sheet of charged impurities with density n 2 
located at ^ = 0, i.e. N[{z) = n 2 S{z). This case was not 
explicitly analyzed by Monroe et al. We therefore numer¬ 
ically integrate Eqns. (4-5) to obtain g for given values 
of 77-2 and n. The resulting density dependence g{n) is 
very weak, with g oc on average. We also find that 
the Dingle ratio rt/rq « 1 for n > 0.9 xlO^^/cm^. 

We now compare the experimental data with these pre¬ 
dictions, starting with the high mobility sample. Eigure 
4 shows g{n) for Wafer No. 16. At low densities, g{n) 
roughly follows the power law expected for remote im¬ 
purity scattering, while for higher densities /i is a weak 
function of n. Superimposed on the data are dashed lines 
showing the expected scaling for remote impurity scatter¬ 
ing and scattering from impurities in the QW. To com¬ 
pare with theory for remote impurity scattering, we set 
Zq = 50 nm which is the SiGe spacer thickness h of this 
sample such that the remote impurities are Al 203 /Si in¬ 
terface charges, as reported by Li et al. [30]. We then ad¬ 
just rii = 3.7 X 10^^ /cm^ to bring theory into agreement 
with the data. Similarly, for scattering from impurities 
in the QW we find reasonable agreement with the data 
when Til = 3.4 xlO^/cm^. Dingle ratio data for Wafer 
No. 16 are plotted in Eig.[^d) and show a linear decrease 
with n over the entire density range. This is broadly con¬ 
sistent with a crossover from remote impurity scattering 
limited transport to local defect scattering-limited trans¬ 
port as n increases. 

In comparison, g{n) is shown for the low mobility sam¬ 
ple (Wafer No. 5) in Eig. 7(a). At T = 0.35 K, the mobil¬ 
ity is weakly dependent on density over the entire density 
range, consistent with scattering from impurities in the 
QW. The dashed line shows the prediction for scattering 
from impurities in the QW taking 77,2 = 1.3 xlO^^/cm^. 
We note that this defect density is 4 times higher than 
the high mobility sample, reminiscent of the factor of 5 
difference between the oxygen contents in the QWs of the 
low and high mobility samples [33|. It is also clear that 
the Dingle ratio is much less sensitive to density, with 
^ 1 “ 2 over the entire density range. The small 
Dingle ratio is consistent with scattering from impurities 
in the QW. 


METAL-TO-INSULATOR TRANSITION 

Eor spin-based quantum information processing, quan¬ 
tum dots are typically operated in the few-electron 
regime El- It is therefore important to characterize 
the strength of the disorder potential in the low electron 
density regime. One important gauge for the degree of 
disorder is the critical electron density, tt^c, for the metal- 
to-insulator transition (MIT) in 2DEGs. Recent experi¬ 
ment and theory have established the MIT as a percola- 









FIG. 8: Wafer No. 16. pxx as a function of T at B = 0 
T forn = 0.34, 0.36, 0.38, 0.40, 0.42, 0.46, 0.51, 0.59, 0.73 
and 0.96 xlO^^/cm^ (from top to bottom). A metal-insulator 
transition occurs at a critical density of ric = 0.46 x 10^^/cm^. 

tion phenomenon, where a fraction of electrons become 
localized by the disorder potential mnn]. As such, Tie 
is an important gauge of the degree of disorder present 
in the QW system, and higher quality samples have been 
demonstrated to display lower values of ric [l3]. Achiev¬ 
ing low Tie is therefore important for few-electron quan¬ 
tum dots, since gate control of electrons may be compro¬ 
mised if the disorder potential is large. For silicon, a MIT 
was first observed in MOSFETs gU 11127], subsequently 
in modulation doped Si/SiGe structures [43l|48l[49] and 
more recently, an undoped Si/SiGe 2 DEG structure [50] . 
and an ambipolar Si-vacuum EET m- In particular, val¬ 
ues of Tie vary greatly in Si/SiGe systems, ranging from 
0.32 to 4.05 xlO^Vcm^ [43l[48H50]. 

The experimental signature for the MIT in 2DEG sys¬ 
tems is a sign reversal in dpjdT, where p is the resistivity 
of the system [44]. Eor n > ric, dp/dT > 0 and the 2DEG 
displays metallic behavior. Eor n < ric, dp/dT < 0 and 
the 2DEG behaves as an insulator. In Eig. we plot pxx 
as a function of temperature for the high mobility sample 
at ten different densities below n = 1.0 x 10 ^^/cm^. We 
observe the following features in this data set: 

1 . At the lowest two densities n = 0.34 x 10 ^^ /cw? and 
0.36 X 10^^ /cm^, dp^^/dT < 0 throughout the measured 
temperature range. In addition, pxx appears to diverge 
exponentially at T < 1 K, indicative of a true insulating 
phase [52] . 

2 . At the next three higher densities n = 0.38, 0.40 and 
0.42 xlO^^/cm^, pxx varies non-monotonically with tem¬ 
perature. While dp^^/dT < 0 at T = 4.2 K, dp^^/dT > 
0 for a small, intermediate temperature range. We note 
that this behavior has also been observed by Lu et al. in 
another undoped Si/SiGe sample [50|, and is common in 
Si MOSEET systems [53] . 



FIG. 9: Wafer No. 16. pxx(B) at T = 0.4 K and n = 6.6x 10“ 
/cm^. The onset of SdH oscillations is observed at Sgff = 0.38 
T. Spin degeneracy is lifted at Bq = 0.88 T and the valley 
degeneracy is lifted at = 1.8 T. 


3. dp^^/dT > 0 at 0.5 K for n > ric, where ric = 
0.46 X 10^^/cm^. 

4. At n > ric, with the exception of n = 0.51 x 10^^ 
/cm^, dp^^/dT > 0 up to a crossover temperature Tc. For 
T > Tc, dp^^/dT < 0. Furthermore, Tc increases with 
increasing n. Das Sarma et al interpreted this behavior 
as a quantum-to-classical crossover naimisiiss]. we 
also find that at T < 1.2 K and n = 0.96 x 10^^/cm^ 
(such that T ^ Tp), Pxx is well approximated by a power 
law relation {pxx — Po)/po = 14.8(T/Tf)+30.7(T/Tf)^/^ 
where po is the zero-temperature resistivity with a value 
of 400 ft based on linear extrapolation. This is also in 
good agreement with theoretical predictions [55] . 

We note that the critical density is comparable to the 
lowest value of 0.32 xlO^^/cm^ that has been reported 
in doped Si/SiGe structures [49], and a factor of 4 lower 
than the value of 1.9 xlO^^/cm^ observed in a previ¬ 
ous work on undoped Si/SiGe structure [50], indicating 
a very low level of disorder in our undoped sample. More 
broadly, the critical density observed in our system lies at 
the lower end of the critical density spectrum [56] . A low 
value of ric = 7.7 xlO^/cm^ was obtained in the GaAs 
system m- 


VALLEY SPLITTING 

Another important figure of merit for the Si/SiGe 
quantum well system is the magnitude of the valley split¬ 
ting. The conduction band of Si has six equivalent min¬ 
ima, or valleys. For Si/Sio. 7 Geo .3 QWs, the in-plane 
strain increases the energies of the four in-plane valleys 
by ^200 meV [T5][36ll38]. The splitting of the two lowest 
lying valleys, Av, is dependent on the magnitude of the 
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vertical electric field in the quantum well and the degree 
of disorder. Degeneracy of these valleys provides an ad¬ 
ditional route for spin decoherence. Measurements of the 
valley splitting in Si MOSFET systems have been exten¬ 
sively performed. Most values range from 0.7 - 1.5 meV 
[SSMT] . with one study reporting a value as large as 23 
meV [62]. In comparison, the valley splitting in Si/SiGe 
systems tends to be smaller, ranging from 0.05-0.3 meV 
[631466] . In this section, we determine Av through careful 
analysis of the SdH oscillations in Wafer No. 16. 

In Fig. we plot with n = 6.6 x 10^^/cm^. 

SdH oscillations are observed above an effective field Heff 
= 0.38 T and have a periodicity of 4 in v. For B > = 

0.88 T, we observe change in periodicity of the SdH oscil¬ 
lations, indicating that spin degeneracy has been lifted. 
The periodicity changes again beyond Hv = 1.8 T, consis¬ 
tent with the lifting of both spin and valley degeneracies. 
We have verified that the spin degeneracy is lifted before 
valley degeneracy using the tilted field method [631465] . 

The energy spectrum of 2D electrons in a perpendicu¬ 
lar field is described by four characteristic energy scales. 
The first is the Zeeman splitting, = gp^B^ where /tb 
is the Bohr magneton and g is the electronic ^-factor. 
The second is E\ = ehB /m* — Fz, which is the Landau 
level spacing minus the Zeeman splitting. The third is 
the valley splitting, Av Finally, clear SdH oscillations 
will only be observed when the Landau level spacing is 
greater than the Landau level broadening T « /i/2rg. 
Spin splitting becomes visible when Ez{Bs) ~ T . Based 
on the effective field at which the SdH oscillations be¬ 
come visible, we estimate T ^ Fi(Heff). We then have 
the relation Ez{Bs) ^ Fi(Heff), allowing us to extract 
g = 3.02. The g-factor is in reasonable agreement with 
the value of ^ = 2.9 ±0.1 at n = 5.9 x 10^^ /cm^ found in 
a previous study M- Based on this experimental value 
of we find T ^ 150 /reV. Finally, the valley degeneracy 
is lifted at the field for which Ay{B = 1.8 T) ^ 150 /reV. 
This value for the valley splitting is substantial and com¬ 
parable to the two-electron singlet-triplet splitting that 
is measured in GaAs quantum dots 00 . 

CONCLUSIONS 

We have measured 26 wafers with different growth 
parameters to identify the dominant mobility limiting 
mechanisms in undoped Si/SiGe QW heterostructures. 
At 4.2 K we find correlations between mobility and oxy¬ 
gen content at the QW as well as the thickness of the 
top SiGe spacer. We have also measured the transport 
lifetime, Tt, and quantum lifetime, Tq, of two Si/SiGe 
QW heterostructures across a wide density range at T 
^ 0.35 K. Based on the density dependencies of the two 
lifetimes, we conclude that the mobility of high quality 
samples with low oxygen content at the QW is mostly 
limited by remote impurity charges. Lower quality sam¬ 


ples with high oxygen content at the QW are limited by 
the impurity charges inside or very close to the QW, con¬ 
sistent with the correlations observed at 4.2 K. To further 
assess the merits of the high mobility heterostructure as 
a platform for spin-based quantum dots, we have mea¬ 
sured a low critical density = 0.46 x 10^^/ cm^ for 
the MIT and a valley splitting Av ^150 /reV. 

While we cannot rule out effects due to other types 
of impurities, our SIMS results suggest that significant 
improvements in the electron mobility may be obtained 
by reducing the level of oxygen content in the Si/SiGe 
heterostructure, particularly in regions close to the QW. 
The SiGe spacer layer thickness can also be increased to 
reduce scattering from charged impurities at the surface 
of the wafer. This second approach has limitations for 
quantum dot devices, as it is desirable to have strong 
in-plane electrostatic confinement, which is harder to ob¬ 
tain in samples with deeper QWs. The magnetic fields 
at which the valley splitting was extracted corresponds 
to a cyclotron radius of ^20 nm, which is a realistic size 
for the lithographic patterning of quantum dots on Si. 
Therefore efforts should also be directed towards reduc¬ 
ing the size of Si quantum dots to emulate high levels of 
magnetic confinement, which yielded large values of val¬ 
ley splitting in this work. Overlapping gate architectures 
may prove helpful to achieve tight electronic confinement 
in the relatively high effective mass Si/SiGe quantum well 
system [67] . 
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